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This article was originally meant to precede a second article [1], where now 
definitions and notations used here are described in more detail. This part will 
here be cut down to a minimum. 

The aim of the investigations is to replace the ordinary treatment of the 
scattering problem by an indirect one. This can be done by providing the 
scattering potential with an adjustable numerical factor which is used as an 
eigenvalue parameter of a discrete eigenvalue problem. 

The term ‘‘resonance method” simply alludes to the structure of the for- 


mulae. 


1, Introduction 


It is well known that in scattering problems an incoming plane wave of 
particles can be separated into waves of different angular momentum, S-, 
P- and P-waves, etc., and that each wave may be treated separately. The 
quantity governing the propOTtiesof the scattered wave is called the asymp- 
totic phase or phase ..Here we shall 

denote it by £,(&) with m l and 

muthal quantum number mmnentum k or energy fc 3 . Smce 

simplicity we are here dealing only with #-waves, 1=0, we may simplify 
the notation into £(&) or and when possible simply f. 

The scattering problem is given by a differential equation 


|^+fc* + ^i7(r)}«/ = 0 


( 1 ) 


* The portion of this work done at the University of Wisconsin Theoretical Chem- 
istry Institute was supported by Grant NsG-275-62{4180) with the National Aero- 
nautics and Space Administration. 

The research carried out at the University of Florida was partially supported by the 
National Science Foundation. 
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with the boundary condition 2/(0) = 0. The function ATJ(r) is the perturbing 
or scattering potential, however of opposite sign, so that — AU(r) is the true 
potential energy. If the solution has the asymptotic form 2/(0“^ sin (fcr + f) 
as r-^oo, | is the asymptotic phase. 

If the potential function does not allow for an elementary solution the 
problem is usually not easily solved directly, even though a number of var- 
iational methods (L. Hulthen, J. Schwinger, etc.) have been developed to 
serve the purpose fairly satisfactorily. 

In this investigation we shall try an indirect method, using the potential 
strength A as a variable eigenvalue parameter, whereas k is so fixed that 
is corresponds to a negative energy state. If k = ix the energy is — x 2 . So 
it is also for k= —ix, but we shall distinguish between the two values so 
that the asymptotic solutions e* XT corresponds to e ±lkr . This means that k- 
values corresponding to bound states are always on the positive imaginary 
&-axis. 

In order not to meet with particular difficulties we shall take U(r) to be 
positive all the way. Then the eigenvalues of A are discrete eigenvalues 
which are usually somewhat easier obtained than the asymptotic phase 
itself by direct means. The calculation of an infinite number of eigenvalues 
A = A n , n= 1,2, is not a serious obstacle because as n -> oo asymptotic 
values are usually fairly easily found. 

The method was first presented at a meeting of the Norwegian Physical 
Society in Trondheim, May 1956, and subsequently published [2] as a 
preprint, which because of several misprints was republished a long time 
afterwards in our Institute Reports [3]. 

As the investigation is presented here it is not much different from the 
institute report, even though former more complete. It was considered to 
use the notations introduced mainly by Jost and Kohn [6] (see references 
in ref. [8]) and adopted in recent works [4], This was finally rejected because 
of some inconveniences and because of the uncertainty as to whether they 
are going to be the final ones [5]. 

2. Definitions and Notations 

We shall consider two particular solutions of eq. (1), + k, r), as our prin- 

cipal solutions. They are defined by their asymptotic values 

u( ± k, r)->e ±ikr as o, (2) 

and hence correspond to the Jost and Kohn functions /( + k, r). They also 
differ from the notation of Ref. 3 where the exponentials are considered 
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split off such that there u( ±k, r)->l. For real k the two functions are con- 
jugate and may be denoted by the shorter symbols u and u*. In the case 
of complex k the shorter notation u may still be used. In particular for k — ix, 
u(k , r) = u(ix , r)->e _KT and hence it is u which has to be used in the eigen- 
value equation 

If we define by u( ± k , r) e ±is — y 2 ± iy l9 

a second set of independent solutions 

y t {r) = \{ue ii +u*e' ii } . 

it is seen that requiring u(k, 0) e l * real means y x (0) = 0. Using the abbrevia- 
tions u( ± k, 0) = u( ± k) we have 


(3) 

(4) 

(4 a) 


u(k) = | u(k) | e w( — ifc) — | w(&) | e* 1 . 


(5) 


u( + k y r) may conveniently be called amplitude functions and u( + k) asymp- 
totic amplitudes. 

The Wronskians of the two sets of functions, defined as W{y, z } — yz — yz\ 
are seen to be 

W{v,yU*} = 2ik, W{y v y^ = k. (6) 

It follows from (4) and (6) that asymptotic values at r = 0 are 


y*(0) = |«(i)|, as r->0 - 


(7) 


3. Solutions by the Theory of Eigenvalues 

Provided U(r) positive in the region 0< oo, eq. (3) has an infinity of 
positive eigenvalues A n {x), n— 1,2, oo. In the case of D r (r)-M7e~ otr as 
r-^oo, e~ ar may be used for a new independent variable transforming the 
fundamental region into a finite region 0 < e~* r < 1, from which it may be con- 
cluded that the eigenvalues increase quadratically with n. 

We now consider u(k) a function of A. Since u(k , 0) = u(k) is bound to be 
zero for any k producing a closed state eigenfunction we must have 


•- 631059 
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«<»*)= n [1- 7-7-:]. 

n-1 L A(*)J 


its value for A = 0 being unity. Provided the A n (x) are known for all real 
values of x and, moreover, may be taken to be analytical functions of x 



u(k) = n 

n = l 

t 1 A n (-ik]’ 

(8 a) 


«(-*)= n 

n = l ! 

t 1 .Mi*)]’ 

(8 b) 

and 


1 -A/A n (ik) 
-A/A n (-iky 

(9) 


4. Zeros and Infinities of the Asymptotic Amplitudes 

If the parameter A fixing the magnitude of the scattering potential is 
given the zeros of u(ix) are found from the equation 

A n (x) = A. (10) 

The A n (x) increase monotonically with x for positive real x. Therefore, if 

A m (0)<A<A m+1 (0) (10 a) 

there will be a positive rvalue for which A m (x) = A, and there will be 
another and higher x - value for which A m _i(«) = A, and so on, the number 
of positive x - values which are zeros of u(ix) becoming ra. In the &-plane 
these zeros k = ix are all on the positive imaginary &-axis. 

In the case of an exponentially decreasing potential of the asymptotic 
form e _ar , x = e~ xr may be used for a new independent variable in a limited 
region 1. If U(r) is positive all the way it is found that asymptoti- 

cally A n (x) is a quadratic function in n. Since A n (x) increases with x it must 
have the asymptotic form 

A„(x) = c(n 4- ax) (n + bx) (10 b) 

with a,b> 0. For larger n therefore u(ix) will have a double set of negative 
zeros corresponding to a double set of k - zeros on the negative imaginary 
axis. In the case of a or b being zero, there is only one such set, the other 
one being skipped away to infinity. The result hence is a finite number of 
zeros of u(k) on the positive imaginary axis, the number depending on the 
magnitude of A , and an infinite number on the negative imaginary A: -axis. 
Similarly the condition for u(k) becoming infinitely large is 
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From (10 b) we see that asymptotically the corresponding x- values are nega- 
tive and the £- values on the negative imaginary axis. Moreover, the infini- 
ties are simple poles except in the case of a — b. 

Similarly, u( — k) has a limited number of zeros on the negative imaginary 
it- axis, all other zeros and poles being on the positive imaginary axis. The 
poles of u( - k) are by some authors called redundant zeros of S( — k) = e~' zi * 
= u(k)/u( — k). 


5. Exemplifications 

Example 1. U(r) - A(A + 1) a 2 (l - tgh 2 or). (12) 

This is a particularly nice example, an Eckhardt potential decreasing 
asymptotically as e' 2ar and with a flat bottom of the symmetric true poten- 
tial — U (r), i.e., U( — r) = U(r). For this reason the transformed equation (3) 

(13) 

with ar^tghar (13 a) 

becomes extremely simple with symmetric and antisymmetric solutions in 
x or r, of which only the antisymmetric ones fulfil the requirement y(0) = 0. 
Hence, using 

A=Z(X - f 1) (13 b) 

for eigenvalue parameter its eigenvalues are 

A n (x) = {2n-\+~^{2n + -^ (n= 1,2,...,°°). (13c) 

The double set of poles are at x m — — mcc and, hence, the poles of u(k) at 

k m =—imcc (m = 1,2, ...,oo). (13 d) 

The double set of zeros of u(k) as obtained from eq. (10) are at x = 
— oc(2n — 1 — A) and x~ — (2n -f A), or 

= — ix(2n— 1 - A), K = - itz(2n -f X) (n == 1,2, °o). (13c) 

If KK3 there is one zero, and if 2m -KK2m+l there are m zeros 
on the positive imaginary fc-axis, all other zeros being on the negative side. 
The asymptotic amplitude u(k) therefore must be 
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u(k) ~ 



(14) 


since the gamma functions of the nominator are accounting for the poles and 
those of the denominator for the zeros, and finally for A — 0, A — 0, its value 
is 1 . By the duplication formula of the gamma function it may also be written 



(14 a) 


In varying A the situation is: At A = 0 the zeros and poles are accurately 
cancelling. As A increases the poles are kept fixed, whereas one set of the 
zeros are creeping up, the other down. For integral A there is a particular 
degeneration, the zeros and poles again cancelling except for some few of 
them on either side, for instance, 


A = 1, 


u(k) — 


k 

k + iot 


A = 3, u(k) = 


k(k — 2 ioc) 
(&^-^a)(fc4-3^a) , 


A = 2, u(k) — 


k — iac 
k + 2ioc 


A-4, 


_ (k ioc) ( k 3 ioc) 

K } (jfc + 2ia)(Jfc + 4ia)’ 


(14 b) 


from which the law for the distribution of zeros and poles in the degen- 
erate case is easily seen. 

The amplitude function u(k , r) is of course easily found in the ordinary 
way as expressed in powers of £(1 — tgha r) or (l-tgh 2 ar). It may be 
given in three useful forms, namely, 


ezp(-fr)»(fc,r)-f'(-J,A + l,l-f.i^2r), (15) 

(ik \ /l — tghar\~ <fc/ * / ik ik 1 — tghar\ 


(15 a) 


exp 


r) u(k, r) 

— tghar\“ ifc/ * / A ik X+l ik ik 2 \ 

-J rri'T- ir«’ 1 -«’ 1 - tgh H (15b > 


■(- 


Putting r = 0 in (15 b) and using the formula 



Resonance Method in Scattering Theory 


9 


F(a , b,c, 1) = 


(c— 1)! (c— 1 — a — 6)! 
(c — 1 — a)! (c— 1 — 6)! 


we obtain the asymptotic amplitude in the form (14 a). 


Example 2. 


U(r) = 


4A 2 aV 2 “ r 

1 — e _2ar ' 


(15 c) 


(16) 


This is a most interesting potential similar to the Yukawa in having a 
pole at the origin. It is also more definite in its results than the above 
Eckhardt potential of Example 1. Putting 


A = X 2 , x = e 2ar , k = ix , 
we have the eigenvalue equation 

J2, 


(16 a) 


f,, d d x 2 1 —x .1 ^ 

| (1 “*>£**- + T -0 ' < 17 > 


Noting the asymptotic factors 

l — e~ 2aT = l —x and e ikr = ^ 


(17 a) 


to be split off before expanding u in a power series it is not hard to see 
that the eigenvalues of A are 

A * {x)= { n + £) “(it) =W ( W + «) (n=1 ’ 2 ’-’ 00) - (1?b) 


Consulting eq. ( 1 1 ), we see that there is now only one set of poles, x = — not , or 
k n = — inot (n = l,2, oo). (17 c) 

According to eq. (10) the zeros are found from 


which gives 



(17 d) 
(17 e) 


Again we see there are only a few zeros on the positive imaginary axis 
giving bound eigenstates. The rest, infinite in number, are on the negative 
side, as are also the poles. There is, however, no general degeneration. Even 
in the case of integral X cancelling of zeros and poles occurs only accidentally, 
even though their asymptotic distributions as oo have the same density. 
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For the purpose of expressing u(k) by means of gamma functions eq. (17 d) 
should rather have been solved with respect to n. We therefore write it in 
the form 


(• + s +, )(* + sr , H <i7f) 

from which we easily obtain the asymptotic amplitude 



Since for A = 0, A = 0, l — ik/2<x, it follows that u(k) = 1. The zeros of u(k) 
are all obtained from the second gamma function of the denominator by 
putting 

ik 

- — (w=l,2, OO), (18a) 

which on squaring l gives (17 d) and (17 e). 

If in this example we go the ordinary way our amplitude function is 
found to be 

u(k,r) = e tkr F^-f x -l, -g + M-f,e-- r ), (18b) 

which, using (15 c), gives u(k) in the form of eq. (18). 

6. Application to a Potential with a Singularity at the Origin 

The singularity of U (r) in Example 2 does not change the origin into a 
singular point of the wave equation. For that purpose a double-pole of 
U(r) is needed, and as an example we shall study the potential function 

U(r) - + l) 1 r + ^(l- tgh 2 r). (19) 

The most interesting feature of this potential is that if we choose 

A = (X + ix)(k + 1 + //) (19a) 

the energy levels of the bound states are the same for all fx and depend 
only on A. The potential is well known from the theory of molecular poten- 
tials [7] under the name of the Poschl-Teller [8] potential, however written 
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in quite a different form. In molecular theory it may be thought of as a 
modification of the Morse potential which gives the same energy levels, 
however a different variation of the rotational constant or mean interatomic 
distance with the vibrational quantum number. This modification is, how- 
ever, for most diatomic molecules the opposite of what is wanted when 
starting with a Morse potential, for which reason the Poschl-Teller poten- 
tial did not prove to be of much use in molecular theory. 

Introducing in the wave equation the new variables 

z = tgh r, k — ix, (l^b) 

it becomes (1 -^£-^^^-^ + ^ = 0. (20) 

This again may be transformed into 

f , d? d 1 — x 2 d 

— (x + p + 1) (x + fi + 2) + ^ j y = 0, (20 a) 

and finally with the new variable 

i7 = l-s*=l-tgh*r (20 b) 

into the hypergeometric equation 








(20 c) 
(20 d) 


The solutions v of this equation are polynomials of degree n— 1, 
n = 1, 2, oo } provided 

A = A n (x) = (2n — 1 + x -f ft) 2 - \ = (2n — 1 + x + ju,) (2n -f x + ft). (20 e) 

Equating this expression with (19 a) we find the ^-values 

x = A + l — 2n, — (A + 2/i-f2 n) (21) 

for which y( 0) — 0. For negative x- values, however, the asymptotic value 
as ig ^(r)— >e |x * r an( j hence must be rejected as a bound state 

function. 
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In order to be able to discuss the scattering properties of our potential 
(19) from the same viewpoint as in Examples 1 and 2 we must redefine the 
amplitude functions u( + k , r) and asymptotic amplitudes. Again we choose 
solution y x (r) and y 2 (r) of eq. (20) with the properties 

yi(r)-*sin(/b*+£), y 2 (r)-> cos (kr + £) (r-*oo), (22) 

with the boundary condition ?/ 1 (0) = 0. Owing to the singularity of the 
endpoint r — 0 of the domain of integration this means that the two functions 
have the asymptotic forms 

yiW-^c 1 r" +1 , y 2 (r)^-c 2 (r-*0). (22 a) 

As before, the amplitude function should have the asymptotic values 
u( ± k } r)^e ±ikr ; however, in order that the asymptotic amplitudes be finite, 
we must now choose 


or 


u(k, r) = tgh^ r(y 2 + iy x ) 

Vi( r ) = tgh” * r • ^ [ ue * - u* e~*] , 
y 2 (r) = tgh~^ r * \ [ue* + u* e~ H ], 


(22 b) 
(22 c) 


The condition ^(0) = 0 now determines the phase shift formally as before 


2 f£ __ ^ 

u(k) 


(22 d) 


however, for ^4 = 0 we have no more £ = 0 or u(k) = 1. Therefore formula 
(8 a) is fallible, which means that, although the zeros of u(k) can be found 
from (10), the infinities of u(k) are not given by (11). The poles of u(k ), as 
might erroneously be concluded from (20 e), are not k — — i(/j, + n) but as in 
Examples 1 and 2 still 

k= —in (n= 1 , 2, ..., oo). (22 e) 


This requires an explanation. Assume the potential U(r) to have the 
asymptotic form e _2ar or for simplicity e _2r as above and take e” 2ar or rather 
x = e~ 2r as a new independent variable. Then for r= oo or x = 0 the wave 
equation has a regular singularity with initial forms of the solution 


x ±*/2 = e =F xr = e ±ikr 


(23) 


Therefore, as x approaches an integer n, the solution beginning with aT* /2 
will have x n “* /2 and higher terms with coefficients of the order of magni- 
tude l/{x — n ), i.e., the solution has poles in x = n, and this means that u(k) 
has poles in k= — in as stated in (22 e). 



Resonance Method in Scattering Theory 


13 


The former general formula (8 a) should therefore be replaced by 

AM -A 


u(ix) = n 


n=i (2 n — 1 + x) (2n~\- x) 


C(n ), 


(23 a) 


where C(n), if needed, is some converging factor independent of x. Intro- 
ducing A and A n (x) from (19 a) and (20 c) the result is 


u(ix) = n 

n -1 


lx A \ lx A+ 1 \ 


where C(n ) has been chosen so that according to the Gauss product formula 
the u(ix) can be expressed by gamma functions, namely, 


u{ix) = 


IMN' 


lx A \ t lx A 4- 1\ 


(23 c) 


From k = ix, x = — ik, and using the gamma function duplication formula, 

2 *( -»)!(-*)! 

(24) 


u(h) = 




This result may be checked by the ordinary method of obtaining u(k , r) 
and u(k), using the differential equation (20 c) in which p 4- 1 is replaced 


by 




\ , d 2 

r / 3 \ 

l d i r , , / „ i\ 

n 


|x+1-|k + - — /iji 

Udi l ~l[ {x+h ~ fl) ~{ A+ l) 

Jr =0 ’ 


(25) 


«/ = tgh"'‘r| 

(1-tgh 2 r\* 12 
i 4 / V ~ 

(25 a) 

Introducing A 

from (19 a) it is found that 



E Ax A x 

C_J (2 _ 2“^2 + 


(25 b) 

Since 

j'l-tgh*^’" 2 

-c + rT‘-" 

(25 c) 


it means that y tgh ;i r has the correct asymptotic form and that we may 
write 
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u{k, r) 

-irP-V'i 


ik _ A _ ik A + 1 

2 2^’ '“¥ + ’T _ 


-ik+ 1,1 — tgh 2 r 


(25 d) 


with the result (24) using (15 c) for F(a , 6, c, 1). 

On comparing (24) with (14 a) for a — 1 we see that the bound states, as 
determined by the second gamma function of the denominator, are exactly 
the same, whereas there is a difference in the asymptotic phase shift amount- 
ing to 

A£ = ar g (f + 1 ! - arg (j + \ + A*) ! • (26) 


For integral /x this means a decrease of 

-arctg^-arctg ... - arctg^, (26a) 

which vanishes as &-> 0. As k-+oo } however, it approaches the value 
Af = — similar to the asymptotic values 


i/ ; 


jlT 

J Ju +i (r )^ sin 



(26 b) 


It is of some interest to compare potentials of the type (19) themselves 
for different ^-values. Beginning with [x = 0 the maximum depth is at r = 0 
and is ^.(A+l). If /x is not too small we make the convenient, however 


insignificant, change of writing 

U(r) = [ - + (A + (1 + i) 2 ] (1 - tgh 2 r). (27) 

The maximum depth is then found to be 

U{rJ-P (27 a) 

at the distance r 0 from the origin given by 

teh ’ r ‘-d^h' ie ^"- (2,b) 

Using 1 - tgh 2 r =» - j, - 4<T 2r - &T 4r + . . . , (27 c) 


1 — tgh 2 r 
tgh 2 r 


4e 


-2 r 


(1 


- A ~ 2r \2 


=’4e-°- T + 8e~ ir + 


(27 d) 
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we find to the relative order of magnitude A 2 /(A + fi 4- £) 2 , 

U (r) = 2A 2 e _2<r ' r,) — A* e~ 4(r-r,) — (e _2<r_r,) - c _4<r_r,) ), (27 e) 

i.e. the potential is approaching the Morse potential. 

If in a similar manner we try to generalize the potential of Example 2 
by an additional term with a double pole at r = 0 no coinciding sets of bound 
states are obtained. 

7. Calculation of the Potential U(r) from the Phase Shift 
for Infinitesimal Perturbations 

If in eq. (9) we take A to be small, then to the first order in A 


or 


if for brevity we write 


t = - y j l 1 1 

g 2irZi\A n (-ik) A n (ih)Y 
i_i y L_1 

2t n-i |^4 n ( — ik) A n (ik) ]’ 

i-( 8 A 

s Wh-«' 


(28) 
(28 a) 
(28 b) 


On the other hand, writing the wave equation of the former y x (r) in the form 
{^ + F + JC7(r)J^(r)=0 (29) 


and combining it with the equation for sin lcr , we find 


Coo 


1c sin £ = A I ?7(r) y x (r) sin hr dr. 

Jo 

As A-M), y(r)->sin hr and sin £->£. Hence 

| = 7 f U(r) sin 2 hr dr— f F(r) 8m2krdr, 

«Jo Jo 

and this equation can readily be converted into 

F(r) = - | °° | sin 2krdk, U(r) = - F'(r), F(r) = f " U(r) dr. 
nj 0 Jr 


Take now the potential U (r) = 


4e 


— 2r 


1 — e 


-2r 


(29 a) 

(29 b) 

(29 c) 
(30) 
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of Example 2, eq. (16), putting for simplicity a = l, we have from (17 b) 
A n {x) = n(n + x) and from (28 a) 


1=^2 


2i n=i [ n(n — i&) n(n + i&) J ' 


(30 a) 


From (29 c), transforming half of the integral into the region — °o ^ k < 0, 
we obtain 

(3o »» 

Closing the path of integration in the upper half plane and contracting it 
into circles around the poles h = in of the second term, the result is 


m- 2 2 ;e- 2 "' = 21og-Ar r , 

71 — 1 tv 1 € 

4e~ 2r 

- F '(r)=U(r) = Y -^- r . 

Similarly, in the case of Example 1 with a = 1 and 

U(r) = 1 - tgh 2 r, A n (x) — (2 n- 1 + x) (2n + «), 

1 


^(r) 


1 00 /*00 1 

= -- 1 \ A 

7T n = l J — oo 1 


ifc) (2w — i&) (2w — 1 + i&) (2n + ifc) 


For the latter term in the bracket we write 

1 


[i-(2ra-l)i] [it -2m] 

and, on proceeding as above, we obtain 


1 [ 1 1 

i [A: — (2 n — 1 ) i k— 2 ni_ 


2e~ 


(30 c) 
(30 d) 

(31) 

e 2ikr dk 

(31a) 

(31b) 

(31c) 
(31 d) 


*(r )- 2 2 (-l)- 1 *-*"-.— zr,-l-tghr, 

n = 1 1 + e 

— F'(r) = 1 — tgh 2 r = U(r). 

8. Derivation of the Coulomb Phase Shift by Confluence 

In the case of a pure Coulomb potential it is well known that the solutions 
of the wave equation are confluent hypergeometric functions, which may 
fairly easily be studied directly. If we write 


( 32 ) 
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for the standard potential function, using ^4£7(r) in the wave equation, it 
is seen, however, that it can be obtained from 


by putting 


A)? a 2 e~ 2ar 4Axe~ 2xr 

(32 a) 

AU (r)- 1 _ e _ 2ar 1 _ e -2*r 

A 

A 2 = — and making a->0. 
a 

(32 b) 


In this way we are able to use all results from the former treatment of the 
potential in Example 2. 

The apparently hardest point to settle in this way is the occurrence of 
the r-dependent part of the phase shift (A/k)\og 2kr. But this is in the 
Coulomb case obtained as an asymptotic part of the solution of the equation 

f d* 2A] _ ikr n x 

< — + 2ik — 4 \ue =0 (32c) 

[dr 2 dr r J 

using only the two last terms, which gives ne~ ikT ~ i(A/k) log r. 

The amplitude function u(k, r) may therefore be defined by its asymptotic 
value 

u(k , r)^ kr+Aik loe2kr) (r->oo), (33) 

and the asymptotic amplitudes as before by u( + k) = u(± k, 0). Then writing 

y i(r) = — { u(k , r) e i$ -u(-k, r) e _< *} , (33 a) 


y x (r)^sm(kr+ (A/k)\og2kr-\- £) as r->° o. (33 b) 


Requiring 2/(0) = 0 we obtain as before the r- independent part of the phase 
shift as 


r 2 ig _ «(-*) 

u(k) ’ 



n(-k) 
u(k ) * 


(33 c) 


It should be noted that the argument kr in the logarithm is suggested in a 
natural way by introducing the variable kr , which also leads to A/k as a 
natural constant of the equation. The factor 2 in the logarithm must, how- 
ever, be justified separately. 

Using the equation 


d 2 „.d 4Aoc e“ 2orr 

dr 2 + tk dr + l-e~ 2 * T 


ue 


,-ikr / 


( 34 ) 


in a similar manner as (32 c) we may argue that the asymptotic factor 
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'2k{l — e -2 “ r )"| u,fc 
2a 


-+(2kr) 


lA/k 


as a-> 0, 


(34 a) 


as obtained from the two last terms of eq. (34) should be taken into the 
asymptotic expression for u(k, r). 

As long as we are dealing with finite a the factor 1 — e -2ar is 1 as r~> oo 9 
so it means that we have multiplied our former u(k , r), and hence also u(k), 
by ( lc/oc ) tAlk . We therefore have to change expression (18) into 



(34 b) 


As a->0 we find 


ik + A 
2a ' ik 


and 



(34 c) 


(34 d) 


From Stirling’s formulae for asymptotic values of gamma functions it is 
found that for large n , nl/(n + £)l = n~ e . Hence 


and 



(35) 


(35 a) 


9. The Coulomb Scattering by the Resonance Method 
The Coulomb potential U(r) = 2/r has an infinity of bound states, 

En = k\= -x\= (n= 1,2, ...,oo). (36) 

71 


Take k = ix to be given and A an eigenvalue parameter, the equation 



% 
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has the eigenvalues A n = nx (n = 1 , 2, . . . , oo ) (37 a) 

and the eigenfunctions 

with the asymptotic values 

y n (0) = 0, yn( r )—>e * T ( ~ 2xr) n . 

Replacing now x by ik and n by A/x^ iA/k we define u(k , r) to be the 
corresponding solution with the asymptotic value 

U (k, r ^ e HkfHAlk) log 2*0 

dropping the factor 

As usual we define two real solutions of the wave equation 


(37 b) 

(37 c) 
the 

(38) 
(38 a) 


Id* 

\dr 2 


d? 2A] 

+ k 2 4 ? y — 0 


(39) 


by Vi = 2i r ) u( — k, r) e l( } — >sin (ir + ^ log 2 hr + (39 a) 

\ {u(lc, r ) e i£ + u( — k, r) e -i£ }-*cos (hr + ~ log 2kr + c j (39 b) 

and obtain from the boundary condition ^(0) = 0 

u(i) = u(l,0). (39 c) 


e a{ = «( - fc) 




The bound states require u(k) — 0 for some positive imaginary k. It follows 
that u(k) must be expressible as an infinite product containing the factors 
1 — A/A n ( — ik). 

It is clear also that u(k) has no poles, contrary to the case of an expo- 
nential potential. The poles of u(k) in eq. (34), for instance, which are 
k= — inoL , are creeping upwards in direction of the origin as a->0. More- 
over, according to (17 e) and (32 b) they are killed by the half-set of zeros, 
n — (A/ol ) * + m, k — — i(Aoc) * — imx, which as a-* 0 have the same distribution 
k — — ima. The other half-set of zeros, n< A/a*, are all of them above the 
origin. As a->0 they become infinite in number and are distributed as 


k n = iA / n. 


(39 d) 
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Since u{k) has only zeros and no infinities we need not worry about the 
converging factors in the above-mentioned infinite product but may write 
at once 


u(k) = 



(40) 


since from 


— = -n (to =1,2, °°), 

ik 


(40 a) 


we obtain all zeros of (39 d). Unknown is still the factor C(A/k), which from 
the transformation properties of the wave equation, is seen to be a function 
of A/k. 

Now the absolute value | u(k) | is easily found from 

kv 

2/ 2 (0) = |m(A:)|, t/(r-*0) = |^y-|, ( 40b ) 


which follows as before from (39 b) and the unaltered value k of the Wron- 
skian. On the other hand, the asymptotic value of (37 b) as r-^0 is 


y n (r->0) = 2 ikr * n\ — 2ikr • 



(40 c) 


and this has to be divided by 2 i from (39 a), by exp ( — \nA/k ) from (38 a), 
and finally multiplied by 


is _ «( ~ &) 

i «(*) i 



(40 d) 


from eq. (40). Hence 


y{r-> 0) = kr • exp 



(40 e) 


which gives from (40 b) 


I u(k) | = 


exp ( — \nA / k) 



(40 f) 


Hence C(A/k) = exp ( - \nA/k) and (40) and (35) are in accordance also 
with respect to the numerical factor. 
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The phase shift is independent of this real factor and could have been 
given from (40) without additional considerations as 

£ = - ar 8 (y) ! = -In log (41) 
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